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Almost all studies of vortex states in helium II have been concerned with either ordered vortex 
arrays or disordered vortex tangles. This work studies numerically what happens in the presence of 
both rotation (which induces order) and thermal counterflow (which induces disorder). We find a 
new statistically steady state in which the vortex tangle is polarized along the rotational axis. Our 
results are used to interpret an instability which was discovered experimentally by Swanson et al. 
years ago and the vortex state beyond the instability that has been unexplained until now. 
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Most configurations of quantized vortices [1] which 
have been investigated in helium II can be grouped into 
two types: ordered vortex arrays and disordered vor- 
tex tangles. An ordered array is made when helium II 
is rotated with angular velocity exceeding a certain 
small critical value. The resulting quantized vortices 
are aligned along the rotation axis and form an array 
whose areal number density is given by Feynman's rule 
Lrot = 2^1/ K, where k = 9.97 x 10~'*cm^/sec is the quan- 
tum of circulation. A spatially disordered tangle is ob- 
tained [2,3] when helium is made turbulent under ther- 
mal counterflow velocity Vns faster than some small crit- 
ical value, or, more in general [4], using grids [5] or pro- 
pellers [6]. In this work we concentrate on counterflow 
because this form of turbulence is simpler (no need to 
worry about large scale motion and effects induced by 
normal fluid eddies). In turbulent counterflow the vor- 
tex line density (length of vortex line per unit volume) 
is Lfiow = iji^ns where is a temperature dependent 
coefficient [3]. This vortex system is almost isotropic, 
provided that one neglects a small anisotropy induced by 
the imposed counterflow [7]. 

An important question naturally arises: what hap- 
pens if vortices are created by both rotation and coun- 
terflow along the rotational axis? We are aware of only 
one experiment which addressed this issue, which was 
performed by Swanson, Barenghi and Donnelly [8] years 
ago. The counterflow channel was mounted on a rotat- 
ing cryostat, so it was possible to create vortices by in- 
dependent combination of rotation and counterflow. The 
absolute value of the vortex line density L was deter- 
mined from the measurement of the attenuation of a 
second sound along the channel, which was calibrated 
against the known density in rotation [9] . Swanson et al. 
found that at slow rotation the critical counterflow veloc- 
ity above which the flow became turbulent was greatly 
reduced. Furthermore two critical velocities Vd and Vc2 
were observed. For Vns < Vd the measured vortex line 
density L = L^^t was independent of the small values 



of Vns and in agreement with Feynman's rule, which 
was evidence for an ordered vortex array. The value 
of Vd was consistent with the critical velocity of a vor- 
tex wave instability which had been first observed by 
Cheng et al. [10] and then explained by Glaberson et 
al. [11]. This instability, hereafter referred to as the 
Donnelly-Glaberson(DG) instability, takes the form of 
Kelvin waves (helical displacements of the vortex cores) 
which are destabilized by the component of the counter- 
flow velocity in the direction along the vortices. 

Unfortunately, it was not clear to Swanson et al. [8] 
what was the nature of the flow beyond the DG insta- 
bility (Vns > Vd)- Their experiment showed that rota- 
tion added fewer than the expected 2Q/k vortex lines 
to those made by the counterflow. Lacking any direct 
flow visualization the nature of the flow was a mystery. 
Theoretically, the linear stability analysis of Glaberson 
determines only the value Vd of the instability of the 
vortex array. What happens beyond the instability can 
be studied only by nonlinear analysis. Thus the pioneer- 
ing experiment of Swanson et al. has lacked a theoretical 
interpretation. 

Motivated by their work, we study numerically quan- 
tized vortices under both rotation and counterflow using 
the vortex filament model. First we prepare an initial 
vortex array in rotation, then we apply a counterflow Kis 
along the rotational axis. We find that when Vns > Vd 
Kelvin waves are excited and grow (DG instability) , and 
the vortex array becomes unstable. Then we show the 
first numerical evidence for a polarized vortex tangle (po- 
larized superfluid turbulence). 

Our work is relevant to other contexts of current in- 
terest, such as superfluid "'He and atomic Bose- Einstein 
condensates [12,13]. In general, our problem is concerned 
with the competition between order (here induced by ro- 
tation) and disorder (here induced by the heat flow). 

For superfluid '*IIe, the vortex filament model is very 
useful, because the vortex core radius aq ^ 10~^ cm 
is microscopic and the circulation k = 9.97 x 10^'* 



1 



cm^/scc is fixed. Neglecting the normal fluid, according 
to Helmholtz's theorem, a superfluid vortex at a point 
moves with velocity which depends on the shape of the 
vortex and on the velocity field Sq at that point which 
is induced by other vortices [14]. Therefore, in order to 
study a rotating tangle, we need to formulate the laws of 
vortex dynamics in a rotating frame [15]. 

The natural way to perform the calculation in a ro- 
tating frame would require to consider a cylindrical con- 
tainer. We do not follow this approach for two reasons. 
Firstly, our formulation is implemented using the full 
Biot - Savart law, not the localized-induction approxi- 
mation often used in the literature. This would require 
to place image vortices beyond the solid boundary to im- 
pose the condition of no flow across it. This has been 
done in cartesian (cubic) geometry, but it is difhcult to 
do in cylindrical geometry. Secondly, the original experi- 
ment by Swanson et al. [8] was carried out in a rotating 
channel with a square cross section. 

In a rotating vessel the equation of motion of vortices 
is modified by two effects. The first effect is the force 
acting upon the vortex due to the rotation. According to 
Helmholtz's theorem, the generalized force acting upon 
the vortex is balanced by the Magnus force as 



PsK{s X So) = -j^, 



(1) 



where F' = F ~ Q ■ M is the free energy of a system in a 
frame rotating around a fixed axis with angular velocity 
n and angular momentum M. Taking the vector prod- 
uct with s', we obtain the velocity Sq. The first term 
F due to the kinetic energy of vortices gives that Biot - 
Savart law, and the second term SI ■ M leads to the veloc- 
ity Srot of the vortex caused by the rotation. The second 
effect is the superflow induced by the rotating vessel. For 
a perfect fluid we know the analytical solution of the ve- 
locity fcub inside a cube of size D rotating with angular 
velocity SI = ^Iz [16]. By taking into account both ef- 
fects, we eventually obtain the velocity sq in a rotating 
frame which is 
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-V Srot + Vcub- 



(2) 



where the first term at the right hand side is the usual 
de-singularization of the Biot - Savart integral which is 
well known in the literature. Here the vortex filament 
is represented by the parametric form s = s{^,t). The 
symbols and Z_ are the lengths of the two adjacent 
line elements after discretization, and the prime denotes 
differentiation with respect to the arc length ^. The sec- 
ond term represents the nonlocal field by carrying out 
the integral along the rest of the filament on which si 
refers to a point. 



At a finite temperature, a quantized vortex is also 
affected by the mutual friction, as the vortex core is 
dragged by the normal flow. Taking this effect into ac- 
count the velocity of a vortex line at the point s is given 

by 

S = So + as' X {Vns - So) - a's' X [s' X {Vns - So)], (3) 

where a and a' are the temperature-dependent friction 
coefficients. 

Two quantities are useful for characterizing the ro- 
tating tangle. The first is the vortex line density L = 
(1/A) J d^, where the integral is taken along all vor- 
tices in the sample volume A. The second is the po- 
larization of the tangle, which we define by < s'^ >= 
{1/AL) J d^s'(^) • z where z is the unit vector along the 
z direction. Note that < .s^ > is unity for a vortex array 
and zero for a randomly oriented tangle. 

We now describe briefiy how to perform the numerical 
calculation [17]. A vortex filament is represented by a 
single string of points at a distance apart. When two 
vortices approach within A^, they are assumed to recon- 
nect. The computational sample is taken to be a cube of 
size D = 1.0 cm. We adopt periodic boimdary conditions 
along the rotating axis and rigid boundary conditions at 
the side walls. All calculations are made under the full 
Biot - Savart law, placing image vortices beyond the solid 
boundaries. Our space resolution is A^ = 1.83 x 10~^ cm 
and the time resolution is At = 4.0 x 10~^ sec. The coun- 
terfiow Vns is applied along the z axis. The normal fluid 
is assumed to be at rest in the rotating frame, and, to 
make comparison with the experiment [8], we use a = 0.1 
and a' = for the temperature T = 1.6 K. 

Swanson et al. [8] found that the first critical velocity 
Vci is proportional to fl^/^, in agreement with the DG 
instability of Kelvin waves. Glaberson et al. [11] mod- 
elled a vortex array rotating at angular velocity as a 
continuum and found that, in the absence of the friction, 
the angular velocity of the Kelvin wave of wavenumber k 
is given by the dispersion relation lv = 2^1 + i/k^, where 
v = (K;/47r)ln(6/ao) and b « L~^/^ is the average dis- 
tance between vortices. This dispersion law has a critical 
velocity Vdg = {^/k)min = 2(2f2i^)^/^ at the critical 
wavenumber koG = -\/20/z/. Therefore, if the axial flow 
Vns exceeds Vdg for some value of k, Kelvin waves with 
that wavenumber k become unstable and grow exponen- 
tially in time. Physically, the phase velocity of the mode 
k is equal to the axial flow, so energy is fed into the 
Kelvin wave by the flow. 

Figure 1 illustrates the DG instability. The compu- 
tation was performed with angular velocity Q = 9.97 x 
10~'^rad/sec, for which Vdg = O.OlOcm/sec. The vor- 
tex lines remain stable when Vns = 0.008cm/sec < Vug 
(Fig. la), while at Vns = 0.015cm/sec > Vdg Kelvin 
waves become unstable and grow (Fig. lb), as predicted. 
Figure Ic shows that Kelvin waves of larger wavenumber 
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become unstable at higher counterflow velocity. This is 
the first numerical confirmation of the DG instability in 
a rotating vortex system. 




(a) (b) (c) 

FIG. 1. Numerical simulations of the Donnelly-Glabcrson 
instability at fl = 9.97 x lO^^rad/sec, T = 1.6K. Each snap- 
shot of the vortex configuration is for the counterflow velocity 
Vns = 0.008cm/sec (a), 0.015cm/sec (b) and O.Ooem/sec (c). 

The key question is: what happens to the vortices be- 
yond the DG instability? Because of the computational 
cost of the Biot - Savart law, it is not practically possi- 
ble to compute vortex tangles with densities which are as 
high (L = C'(10'*)cm~^) as those achieved in the exper- 
iment. Nevertheless the following numerical simulations 
can shed light into the physical processes involved. The 
time sequence contained in Fig. 2 illustrates the evolu- 
tion of a vortex array at 17 = 4.98 x lO^^rad/scc under 
the counterflow Vns = 0.08cm/sec. Figure 2a shows the 
initial N = 33 parallel vortex lines; they have been seeded 
with small random perturbations to make the simulation 
more realistic. As the evolution proceeds, perturbations 
with high wavenumbers are damped by friction, whereas 
perturbations which are linearly DG-unstable grow expo- 
nentially, hence Kelvin waves become visible (Fig. 2b). 
When the amplitude of the Kelvin waves becomes of the 
order of the average vortex separation, reconnections take 
place (Fig. 2c). The resulting vortex loops disturb the 
initial vortex array, leading to an apparently random vor- 
tex tangle (Fig. 2d). After the initial exponential growth 
(which is predicted by the theory of the DG instabil- 
ity) , nonlinear effects (vortex interactions and vortex re- 
connections) become important and nonlinear saturation 
takes place. The two effects are apparent in Figure 3a, 
which shows the initial exponential growth of L (linear 
instability) and the successive equilibration to a statisti- 
cal steady state (nonlinear saturation). 

Looking carefully at the saturated tangle in Fig. 2d 
wc notice that there are more loops oriented vertically 
than horizontally. Figure 3b shows how the polarization 

< s'^ > changes with time in the calculation presented 
in Fig. 2. During the exponential growth (linear phase) 

< s'^ > decreases from unity (ordered vortex array), but 
it never becomes zero (random tangle), settling instead 
to the finite value < s'^ >« 0.5 as soon as nonlinear 
saturation takes place at about t ~ lOOsec. 




(a) (b) 




(c) (d) 
FIG. 2. Numerical simulation of rotating turbulence at 
T = 1.6K, n = 4.98 X 10"^rad/sec and Vns = 0.08cm/sec. 

Computed vortex tangle at the following times: (a): t=Osec; 
(b): t=12sec; (c): t=28sec; (d): t=160sec. 




t (sec) t (sec) 

(a) (b) 
FIG. 3. Time evolution of L (a) and < s'^ > (b) at 
T = 1.6K and V„s = 0.08cm/sec for fi = 4.98 x 10"^rad/sec. 

Figure 4 shows the calculated dependence of the vortex 
line density L on the counterflow velocity Vns at different 
rotation rates Q. The dependence of L on Vns is similar 
to what appears in the Figure 1 of the paper by Swan- 
son et al [8]. The only difference is that the scale of the 
axes in the paper by Swanson et al. is bigger - in this 
particular figure they report vortex line densities as high 
as L Ki 2500cm~^, whereas our calculations are limited 
to L w 80cm~^. Despite the lack of overlap between the 
experimental and numerical ranges, there is clear quali- 
tative similarity between the figures. It is apparent that 
the critical velocity beyond which L increases with Vns 
is much reduced by the presence of rotation, which is 
consistent with the experiment. Figure 5 shows the cal- 
culated polarization < s'^ > as a function of counterflow 
velocity Vns at different rotation rates It is appar- 
ent that the polarization decreases with the counterflow 
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velocity and increases with the rotation. 
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FIG. 4. Vortex line density L vs V^^ aX T = 1.6K for 
n = (write circle), Q, = 9.97 x 10~*rad/sec (black circle), 
n = 2.99X 10"^rad/sec (triangle) and Q = 4.98x 10"^rad/sec 
(square) . 
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FIG. 5. Tangle's polarization < s'^ > vs Vn, at T = 1.6K 
for n = 9.97 X lO^^rad/sec (circle), O = 2.99 x 10"^rad/sec 
(triangle) and n = 4.98 x 10~^rad/sec (square). 



In conclusion, we have studied numerically vortex tan- 
gles under the effect of rotation for the first time. At 
velocities higher than the onset of the DG instability we 
have determined the existence of a new state of superfluid 
turbulence (polarized turbulence) which is characterized 
by two statistically steady state properties, the vortex 
line density and the degree of polarization. Although the 
computed range of vortex line densities does not overlap 
with the much higher values obtained in the experiments, 
we find the same qualitative dependence of vortex line 
density versus counterflow velocity at different rotations. 

Further work will investigate other aspects of the prob- 
lem, particularly the nature of Vc2 and what happens 
at very high counterflow velocities. We also hope that 
this work will stimulate more expc;rinic;nts. For example, 
it should be possible to observe the polarization of tur- 
bulence by using simultaneous measurements of second 
sound attenuation along and across the rotation axis. 
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